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We discuss how to construct stable multidimensional extensions of one-dimensional dark solitons, 
the so-called soliton stripes, in two-species Bose-Einstein condensates in the immiscible regime. We 
show how using a second component to fill the core a dark soliton stripe leads to reduced instabilities 
while propagating in homogeneous media. We also discuss how in the presence of a trap arbitrarily 
long-lived dark soliton stripes can be constructed by increasing the filling of the dark stripe core. 
Numerical evidences of the robustness of the dark soliton stripes in collision scenarios are also 
provided. 
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I. INTRODUCTION 

One of the most remarkable achievements in quantum 
physics in the last decade was the Bose-Einstein con- 
densation of ultracold alkaline atoms [l[ . Because of the 
interatomic interactions present in Bose-Einstein conden- 
sates (BECs), they have a high potential for supporting 
quantum nonlinear coherent excitations of which many 
types have been already been experimentally observed: 
dark @, 0|, bright 3 and gap solitons @, vortices @, 
vortex ring s Q , shock waves [8| , vector solitons Q , and 
others [lfj. Many other types of nonlinear excitations 
have been theoretically predicted to exist (see e.g. Refs. 
[Til [l2j and references therein) . 

Although solitons have attracted a lot of attention be- 
cause of their robustness, they do not survive in effective 
dimensions higher than one: bright solitons are unsta- 
ble to blow-up [TH, [3] and dark soliton stripes decay 
into point vortices due to the so-called snake instability 
• The development of this instability in the context of 
BECs has received a lot of attention [15J . Other multidi- 



mensional extensions of the dark soliton, such as the ring 
dark soliton , are also unstable [TtJ • 

In niulticomponent Bose-Einstein condensates the ex- 
tra degrees of freedom provided by the existence of an- 
other component allows for the existence of many differ- 
ent types of solitons. In the context of BEC there has 
been a strong interest on direct extensions of the one- 
dimcnsional (ID) solitons, such as dark-dark, or bright- 
bright solitons [l§ - l22j and also on many other types of 
nonlinear waves that do not have direct ID analogues 
such as dark- bright solitons (23rl25j . soliton molecules 
domain walls (UE^, symbiotic solitons USUI, 
solitons on plane wave solutions [32 |. and supersolitons 
[33l | . However, as it happens in the scalar case most of 
these solutions arc unstable in multidimensional scenar- 
ios. 



In this paper we will study how adding a second com- 
ponent allows for the construction of stable multidimen- 
sional extensions of the dark soliton stripe thus leading 
to stable dark-bright soliton stripes (DBSS). The idea, 
as in the dark-bright soliton [23H25| concept, is based 
on using a second immiscible condensate to fill the dark 
stripe thus increasing the rigidity of the coupled system. 
Although different types of nonlinearly coupled struc- 
tures have been studied in niulticomponent BECs (see 
e.g. Refs. [U,[3^] and references therein), to our knowl- 
edge this possibility has not been studied in the context 
of niulticomponent trapped BECs. 

The structure of the paper is as follows. First, in Sec. 
HT1 we present our model problem and the basic facts in 
one-dimensional scenarios. In Sec. [TTT] we discuss the 
behavior of dark-bright soliton stripes in homogeneous 
media and discuss their linear stability. In Sec. HVI wc dis- 
cuss the novel features due to the presence of transverse 
trapping and show how the combination of the stripe 
filling due to a second component plus the transverse 
trapping leads to extremely robust dark-bright soliton 
stripes. This robustness manifests also in dark-bright 
soliton stripe collisions, studied in Sec. [V] Finally, in 
Sec. IVII we summarize our conclusions. 



II. MODEL EQUATIONS AND DARK-BRIGHT 
SOLITONS 

A. Model equations and parameters 

We will consider a quasi-two dimcnsinal dilute binary 
mixture of BECs in a pancake trap described in the mean 
field limit by the adimcnsional equations 



dt 



~A + V{x,y) 



£j&fcW a nfc. (i) 
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where j — 1,2 and A = d 2 j dx 2 + d 2 / dy 2 . The wavefunc- 
tion ipj is related to the one in physical units (denoted 
by tilde) in the full space through 
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where \Fj are normalized to the numbers of atoms Nj: 
f\*j\ 2 dc = Nj, and ip (z) = (1/V) 1 / 4 cxp(-z 2 /2), r = 
r/Z z , t = tuj z where l z = \/ h/muj z and lu z is the trap 
frequency along z direction corresponding to a trapping 
potential V(z) = muj 2 z 2 /2 in the original model written 
in physical units. Finally, the nonlinear coefficients are 
defined through gjk = Ojfc/ciii, where an is the intra- 
component scattering length of the first component (as- 
sumed, as all of the others, to be positive). 

For simplicity, all the examples to be discussed in this 
paper will refer to the case of equal inter- and intra- 
species interactions: gn = gi2 = g 2 \ = .922 = 1, that 
is not far from several physically relevant situations. We 
have checked that our ideas to be presented in this paper 
remain qualitatively valid for a wide range of parameter 
combinations corresponding to two-component BECs in 
the immiscible regime. 



B. Dark- bright solitons 

Let us first consider the situation V(x, y) = and take 
initial data of the form 



■4> 2 {x,y,t) 



^d{x 1 t)f{y), 
4'b(x,t)f(y), 



(3a) 
(3b) 



where ipd, 4>b solve the coupled NLS equations [2c 
.dipd ld 2 ifj d 



'St = -2W + 0*" 12 + W - ") *" < 4a > 



Eqs. (Q| have exact dark-bright soliton solutions [2^ 

i/jd(x,t) = y/p[i since + cosatanh/3 [x — X(t))] ,(5a) 
Mx,t) = A 6 sech/3(a;-X(i))e 4W+;r/3tanQ+ni ' t) , (5b) 



where A b = = ^ pcos 2 a + (N b /A) 2 - N b /4 

with N b = J \ip b \ 2 dx and n b = f3 2 (l - tan 2 a)/2-6 (in 
what follows we will take 6 = 0). The soliton position is 
determined by the equation X(t) = xq + 1/3 tana. In the 
scalar case the parameter a is related to the dark soliton 
velocity as v = y^psina, but in the vector case due to 
the presence of the bright component the velocity of the 
vector soliton for the same a is slower, and given by the 
formula Vd b = (3 tan a. 

The bright component also has its own global phase 
cj) that is relevant only when interactions between two 
bright solitons are considered. 



It is worth to mention that this dark-bright soliton 
exists only in the immiscible regime of repulsively inter- 
acting systems and has nothing to do with the bright 
solitons existing in systems with attractive interactions. 
Through this paper when speaking of "bright solitons" or 
"the bright component" we will refer to the repulsively 
interacting component that is hosted by the dark soliton. 



III. DARK-BRIGHT SOLITON STRIPES AND 
THEIR STABILITY 

Taking f(y) = 1 in Eq. ^ an exact solution of Eqs. 
([T]) with appropriate boundary conditions in two spatial 
dimensions is given by Eq. (|5|) . 

In the scalar case, this dark soliton stripe is known to 
be unstable [36j . Here we will check stability of the DBSS 
considering Eqs. (JTJ) and taking the bright component as 
unperturbed, namely ip b (x,y) = ip b o( x ) while the dark 
component is subject to the transverse perturbation 

Mx,v) = 1>do(x) + a{x)e^ k y +mt + b*(x)e ik y~ mt . (6) 




FIG. 1: Surface plot of the critical wavenumber k c for which 
the instability appears as a function of the relevant parame- 
ters (w,Nt). 

Substituting into Eqs. fTJ) in the frame mov- 
ing with the velocity of the vector soliton, Vdb, and de- 
composing a(x) and b{x) into real and imaginary parts 
a = h + ih and b = fx - if 2 we get 



F + LF = 0, 



(7) 



where F 



L = Li 



1 

-1 

iQ~ 



and 

iQ- 



being 



2 w xx 1 



Li = v db ad : 
P 

q ± = 4(A±(4f 



2 \tpdo\ 2 + \ipbo\ 2 - P, 



(8) 



(9a) 
(9b) 
(9c) 
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Rewriting the unperturbed solution given by Eq. ([5]) in 
the moving frame, defining w = y/pcosa, rewriting (3 
and Vdb in terms of w and inserting tpdo and ipbo in Eqs. 
© we find 



L = Li - 



' Qw 2 ~N b f} 
2 cos. 

-2 UW tanh(/3x) zZlh-TA 



^^-2w 2 -2vwiwih.(fix) 



2v 2 



(10) 



The eigenvalue problem can be rewritten in the form 

LF = -\p 2 F 1 (11) 

where \p 2 = iila — \k 2 and p 2 plays the role of a new 
eigenvalue. Finally, after introducing the rescaling 

{n/p 2 ,k//3,/3x, P /p,v//3,w//3} = {n,~k,x,p,v,w} 

(12) 

and looking for rescaled eigcnfuctions of the form 
h=e^ gi {x), f 2 = e i « s g 2 (x), 



(13) 



where g\(x) and g2(x) are unknown functions which are 
constant in the limit \x\ — > oo, of the rescaled eigenvalue 
problem 




FIG. 2: Density plots of the dynamics of a perturbed dark 
soliton stripe with a homogeneous background, (a-c) Scalar 
case, i.e. Nf, = 0. Shown are the densities \tp(x,y,t)\ 2 for (a) 
t = and (b) t = 180 and (c) the phase arg(^) for t = 180. 
(d-f) Vector case with a bright component Nb = 3. Shown are 
the densities \ipd{x, y, t)\ 2 for (d) t = and (e) t = 570 and (f) 
the second component \i/jb(x, y, t)| 2 for t = 570. Parameters 
are p = 1, A = 0, « = and the spatial region shown in 
the all the subplots is (x,y) e [-40,40] x [-40,40]. The 
initial stationary profile was perturbed by taking i^j = %f)j (t = 
0)(1 + 0.005 cos(O.ly)) 



LF = ~fF, 



we obtain, after some algebra 
-2w 2 



-2vw — iqv/w 



—2vw + iqv/w 
2v 2 - \q 2 + \p 2 



(14) 



= 0. 



(15) 



[3fd) (notice the wider dark stripe in comparison with 
Fig. [5{a) due to the presence of the bright component) 
becomes unstable for t ~ 570, these times being three 
times longer than the instability times found for Nf, = 0. 
By increasing the number of atoms in the bright compo- 
nent it is possible to get stable dynamics for much longer 
times. 



Solving the eigenvalue problem (fl"5"|) we get the relation 
between p and q 

p 2 = 2{i) 2 + w 2 ) +q 2 ± 2V(w 2 + w 2 ) 2 + v/(wq 2 ). (16) 

It can be shown (see [HI) that an eigenvalue from the 
discrete spectrum of the operator L corresponding to a 
bound state with Q = is q = i what leads to the follow- 
ing critical wavenumber below which the long- wavelength 
instability appears in the non-scaled variables 



(3 2 - 2 + 2v//3 4 (l - l/«; 2 ) + 1. 



(17) 



In Fig. Q] wc plot the dependence of the critical 
wavenumber k c on the relevant parameters (ui, iVf,). It is 
obvious that the presence of a nonzero bright component, 
although unable to completely inhibit the instability, re- 
duces significantly the range of unstable wavenumbers. 

This analysis is confirmed by direct numerical simu- 
lations of Eqs. ([1]) for the evolution of a black stripe 
on a homogeneous background. Some results are sum- 
marized in Fig. [2j It is clearly seen how introducing a 
bright component leads to a substantial increase in the 
time at which the snake instability sets in. For a small 
bright component with Nf, — 3, the DBSS shown in Fig. 



IV. WAVEGUIDE CONFIGURATIONS 

Let us now consider more realistic scenarios by includ- 
ing a transverse trapping potential in a waveguide geom- 
etry taking V(y) = n 2 y 2 /2. 

Due to the spatial inhomogeneities the analytical sta- 
bility analysis is intractable and we will base our analysis 
on numerical simulations. We will take the transverse 
part of the solution f(y) in Eq. (J3j) as given by the solu- 
tion of the nonlinear spectral problem for (/x, / M ) 



2 dy 2 



1 



2 2 

k y 



f - pf = o, 



(18) 



with boundary conditions /(0) = 1, / y (0) = and 
f(y) = 0. A typical initial profile for the waveg- 
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uide configuration calculated from Eqs. ([5]) and (JT8J) is 
shown in Fig. [3][a). 

We have used these solutions as initial data for a nu- 
merical study of the stability of trapped DBSSs. In all of 
our simulations to be presented in this section, the con- 
densates are far from the ID limit since in that limit dark 
solitons are stable and our dark stripes are all unstable 
when their cores arc empty. In this case, and due to the 
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FIG. 3: Simulation of the evolution of a DBSS. (a) Two- 
dimensional surface plots of a DBSS with initial data, 
\ip d (x,y,0)\ 2 , \ip b (x,y,0)\ 2 given by Eqs. (0, with / given 
by Eq. (|18p for parameter values p = 1, n = 0.4, JVf, = 3, 
cf> = 0, and w = —0.6. (b-i) DBSS propagation for parameter 
values p = 1, k = 0.1, i! = —0.6. Shown are density plots of 
\ip d {x,y,t)\ 2 for (b-d) N b = and times (b) i = 0, (c) t = 50, 
(d) t = 100 and for N b = 3 (e-g) for times (e) t = 0, (f) 
t — 100 and (g) t = 200. Density plots of the bright compo- 
nent in the multicomponent case are also shown in subplots 
(h-i) for t = (h) and t — 200 (i). The spatial region shown 
in subplots (b-i) is x G [-80,80],?/ G [-20,20]. 



previously discussed enhancement of the stability due to 
the filling of the dark soliton stripe it is possible to make 
an unstable scalar configuration to be stable, leading to 
stable multidimensional DBSSs. Fig. G3c-g) shows an 
example of this behavior. 

In Fig. Elb-d) we can see that the scalar dark soliton 
stripe becomes curved when propagating in the waveg- 
uide, a well known effect due to the dependence of the 
soliton velocity on the (now inhomogeneous) background 
amplitude. Thus in the low amplitude regions the stripe 
moves more slowly and transforms into smaller amplitude 
gray solitons accompanied by the formation of an array 
of vortices on the tails of the radial density distribution. 
This behavior is confirmed by the analysis of the phase 
distribution of the dark component. 

It is known that the transverse trapping plays a stabi- 
lizing role in the dynamics of multidimensional solitons 
that in a certain tight-binding regime become effectively 
ID and thus stable [l5[. This fact is also in agreement 
with the ideas of Ref. [23| where it was suggested that 
in higher dimensions dark-bright solitons should be more 



FIG. 4: Density plots of the results of the simulation of the 
transition to more stable behavior of DBSSs as the number of 
particles in the bright component Nb increases for parameter 
values k = 0.2, (a-c) Nb = 0, (d-f) N b = 1, (g-i) N b = 3, (j-1) 
N b = 5. The times are indicated on each subplot. The spatial 
regions shown in all of the subplots is (x,y) G [—20,20] x 
[-10,10]. 

stable than pure dark solitons, even in traps that do not 
attain the quasi-lD regime. However our results go be- 
yond that since we have tested many different config- 
urations including very weakly trapping potentials and 
have found that it is allways possible to find a filling 
that makes the DBSS stable for very large times. The 
transition from the unstable regime to the stable one as 
the number of particles in the bright component is in- 
creased is very clear, the number of vortices generated 
being smaller and the instability setting in later as a func- 
tion of Nb. 

An example of that behavior is shown in Fig. 2J It can 
be seen how for a given trapping, the scalar dark soliton 
stripe becomes unstable and breaks into four point vor- 
tices [Fig. U a-c)], as the bright component is increased 
to N b = 1 [Fig. Hd-f)], N b = 3 [Fig. Hg-i)], and N b = 5 
[Fig. |Hj-l)], the number of vortices arising after the in- 
stability develops are reduced to three, one and none, 
respectively. The times required for the instability to set 
in become larger with increasing Nb as discussed previ- 
ously. For N b = 5 the DBSS is stable for the very long 
times used in our simulations. 



V. INTERACTION BETWEEN DARK-BRIGHT 
SOLITON STRIPES 

Robustness in soliton-solitons interaction is one of the 
key elements of these nonlinear waves. In our case, de- 
spite our 2D dark-bright stripe structures studied are not 
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strict (intcgrablc) solitons, they display collisional prop- 
erties that support their robustness. In this section we 
will consider fast and slow DBSSs collisions as well as 
their interactions in ring trap configurations. 



A. Fast collisions 

First, in fast head-on collisions DBSSs show clearly 
their robustness. We have simulated many collision sce- 
narios and the typical outcome is a pair of outgoing DB- 
SSs with the usual time delay due to the intreraction. A 
typical simulation is shown in Fig. [5] 







■1 


1 (e) 







FIG. 5: Interaction of two DBSSs separated initially in space 
by a distance £ = 40. Parameters are the same for both soli- 
tons: /0 = 1, \a\ = 0.6, k = 0.1, A<£ = 0, £ = 20, N b = 3. 
The upper/bottom row corresponds to the dark/bright com- 
ponents for (a) t = 0, (b) t = 70, and (c) t = 140, respectively. 
The spatial region shown is x £ [—50, 50], y G [—20, 20]. 



J 



B. Ring trap configurations 



As another demonstration of the robustness of the DB- 
SSs, we have simulated numerically their propagation 
and collisions in toroidal traps [13, [H[ . These configura- 
tions have been suggested to provide interesting test beds 
for the propagation of different types of solitons (39M42| . 
In our calculation we adapt the waveguide configuration 
of the trap by using a trapping potential of the form 



V{x,y) = v 2 (r-Rf/2, 



(19) 



where R is ring radius and r = \Jx 2 + y 2 . A typical 
example is shown in Fig. [5] were both the stable propa- 
gation of two dark-bright solitons along the toroidal BEC 
and their stable behavior after the collision, returning to 
the initial state for t = 1000 arc shown. 




£=0 



£=200 



£=500 




£=700 



c > 

£=1000 



(e) 



FIG. 6: Collision of DBSSs in a toroidal trap. The soliton and trap parameters are k — 0.2, R = 40, ATj = 6, v = 0.9. Shown 
are surface plots of the dark (upper row) and bright components (lower row) on the spatial region (a:, y) £ [—50, 50] x [—50, 50] 
for times: (a) t = 00, (b) t = 200, (c) t = 500, (d) t = 700 and (e) t = 1000. 



C. Slow collisions 

As it can be expected a priori, slow (low energy) DB- 
SSs collisions lead to much more complex dynamics. To 
show it in this subsection we will consider the interac- 
tion of two DBSSs with zero initial velocity (black dark 
components) initially separated by a distance £. In this 



case the phase difference between the bright components, 
i.e. A</> = 0i — 02 plays an important role. When the 
phase difference is zero (A<f> = 0) the bright components 
in the DBSSs repel each other while when this difference 
is A(f> = 7r they have an attractive interaction as it hap- 
pens in the scalar case. As to the dark components, their 
interaction is always repulsive. 
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We present an example of our numerical experiments 
in Fig. [7J For small bright components, iV& = 1 (see Fig. 
[7]) and due to the repulsive force between the dark soli- 
ton stripes two pairs of vortices are excited and counter- 
propagate (see the outer part of the figure) while another 
two pairs also repel each other but move more slowly. In- 
creasing the bright component to N}, = 5 (Fig. [3]) leads 
to a more stable propagation although still some profile 
oscillations and vortex generation are observed. Finally, 
using a larger bright component iV& = 8 leads to a per- 
fectly stable propagation of the DBSSs (Fig. [SJ. 






FIG. 7: Interaction of two DBSSs separated initially in space 
by a distance £ = 20. Parameters are the same for both 
solitons: p = 1, a = 0, k = 0.1, A<?!> = 0, £ = 20. N b = 1. The 
spatial region shown in all of the subplots is a; £ [—80, 80], y £ 
[-20,201. 






FIG. 8: Same as in Fig. [7] but for N b = 5. 

By changing the interaction between the bright compo- 
nents from repulsive (A</> = 0) to attractive (A<p = n) it 
can overcome the repulsive interaction between the dark 
components leading to a direct collision between the DB- 
SSs. The outcome of such collisions is shown in Fig. [TO] 




FIG. 9: Same as in Fig. but for N b 



We will study two situations, the first summarized in Fig. 
T0ta-c) where the bright component is below a critical 
value of Nf, and the attractive interactions are not suf- 
ficient to overcome the repulsive force coming from the 
interaction between the dark stripes. The second case 
studied is summarized in Fig. I10f d-f). where the attrac- 
tion leads to head-on collision between both DBSSs. In 
both cases the interaction leads to the destabilization of 
the dark soliton stripes. This is due to the fact that 
in these types of slow collisions marginally stable multi- 
dimensional solitons have a lot of time to interact and 
thus the outcome in many cases is the destabilization 
of the structure [26|, HH . To overcome instability effects 
we have increased the bright component and observed 
periodic collisions between two DBSSs. In Fig. QT] we 
show the first period of the oscillations observed in that 
scenario. However, in this very unstable situation even 
with this strong bright component after the third colli- 
sion (period) the instability sets in and DBSSs break into 
two pairs of vortices propagating outwards (a similar phe- 
nomenon with different types of solitons was discussed in 
Ref. H|). 



VI. CONCLUSIONS 

In this paper we have discussed the stabilization effect 
of a bright component on dark soliton stripes in multi- 
component BECs in the immiscible regime. 

We have shown using analytical tools and numerical 
simulation that dark-bright soliton stripe complexes are 
less unstable in homogeneous systems than their dark 
soliton stripes counterparts without filling. We have also 
studied the combined effect of trapping and filling with a 
second component, that leads to stable dark-bright soli- 
ton stripes. Finally we have presented some examples 
of the robustness of these dark-bright solitonic stripes in 
collision scenarios. 
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FIG. 10: Interaction of two DBSSs separated in space by a 
distance f = 20. Parameters are the same for both solitons: 
p = 1, a = 0, k = 0.1, Acj> = tt, £ = 20. (a-c) 7V 6 = 3.6 
(f = 0; 150; 200), (d-f) N b = 3.7 (t = 0; 220; 360). The spa- 
tial region shown in all of the subplots is x £ [— 50,50], y £ 
[-20,20]. 




.... 



FIG. 11: Same as in Fig. HD1 but for N b = 8. Shown are 
density plots of the dark component for times (a) t = 0, (b) 
t = 430 and (c) t = 600. 
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